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Abstract 

We study the evolution of the Earth collision probability of asteroid 2008 TC3 using 
a short observational arc and small numbers of observations. To assess impact 
probability, we use techniques that rely on the orbital-element probability density 
function characterized using both Markov-chain Monte-Carlo orbital ranging and 
Monte-Carlo ranging. First, we evaluate the orbital uncertainties for the object 
from the night of discovery onwards and examine the collapse of the orbital-element 
distributions in time. Second, we examine the sensitivity of the results to the 
assumed astrometric noise. Each of the orbits obtained from the MCMC ranging 
method is propagated into the future (within chosen time bounds of the expected 
impact), and the collision probability is calculated as a weighted fraction of the 
orbits leading to a collision from the Earth. We compare the results obtained with 
both methods. 

Keywords: Asteroids, dynamics, Collisional physics, Celestial mechanics, Impact 
processes, Near-Earth objects, Orbit determination 



1. Introduction 

The near-Earth object 2008 TC3 was a small fast rotating asteroid (discovered 
by RA. Kowalski, Catalina Sky Survey), which impacted the Earth on October 
7, 2008. This extraordinary object was discovered and then tracked for about 
20 h before its collision. It was approaching at rather small speed of about 7.6 
km/s. Periodic light variation indicated the fast rotation of the object (period 
about 1.5 min) in a complex, non-principal axis ("tumbling") state [18]. 2008 TC3 
entered the atmosphere at 02:46 UTC and exploded 37 km above northern Su- 
dan. The first prediction (S.R. Chesley, JPL) indicated explosion and a rain of 
small meteorites. As expected number of small fragments was found. A total 
of about 600 fragments (10.5 kg) were collected [in]. The meteorites were later 
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identified as ureilites. Collision probability computations and impact prediction 
were performed at the University of Pisa (Italy) and Jet Propulsion Laboratory 
(USA). Preliminary orbit computation immediately indicated collision with the 
Earth, and subsequent "follow-up" observations were performed, finally gathering 
a total of 570 astrometric observations from 26 different observatories. NEODys 
(Pisa) lists five other independent possible previous close approaches of 2008 TC 3 
with the Earth and one with Mars. Four of which have the probability approach- 
ing 100% and one approaching 50%. In Sect. 2 we describe the methods used for 
computing orbits and collision probability. In Sect. 3 we present our results and 
Sect. 4 contains conclusions. 

2. Orbits and collision probability computation 

2.1. Orbital inverse problem 

An outline of the Bayesian treatment of the orbital inverse problem can be 
found, for example, in J2U1 EH HI I2"T] . The solution to the orbital inverse problem 
can be written posteriori probability, using Bayes' theorem: 



• p pr (P) is the a priori probability density function (p.d.f.), 

• p e (A\E'(P)) is the observational error p.d.f. (usually being assumed Gaus- 
sian), evaluated for the observed-minus-computed (O-C) residuals A^(P), 

• C = (J p(P, il))dP)~ l is a normalization constant, 

• P denotes osculating elements (x, y, z, x, y, z) T or (a, e, i, Q, u, M ) T at epoch 



• if) denotes a set of astrometric observations, consisting of R.A., Dec. pairs: 
ijj = (a^di, ...;a N ,5 N ) T . 

To keep p p invariant in transformations from one orbital element set to another, 
one can use, for example, Jeffreys' a priori p.d.f. [9]: 



p p (P) = C Ppr (P)p(ij | P) 



Qv(P)p e (A*(P)), 



(1) 
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where £ _1 is the information matrix (or the inverse covariance matrix) eval- 
uated for the local orbital elements P, <fi contains the partial derivatives of R.A. 
and Dec. with respect to the orbital elements, and A is the covariance matrix for 
the observational errors. For details of asteroid orbital inversion theory, refer to 
the extensive literature available [201 H31 ED ED E] • To characterize the complicated 
a posteriori target densities, we use Monte-Carlo methods. In particular, we focus 
on the Markov-chain Monte-Carlo ranging technique. 

2.1.1. Markov-chain Monte-Carlo ranging 

Markov-chain Monte-Carlo (MCMC) ranging [15J is a nonlinear orbital sam- 
pling technique that combines the Monte-Carlo (MC) ranging method and Metropolis- 
Hastings (MH) algorithm. In contrast to MC ranging, in MCMC ranging the sam- 
pling is guided, with a help of arbitrary proposal density (here we use multi-normal 
proposal density) for inversion parameters. Subsequent orbits result in a chain 
(Markov chain) of samples which are not independent. As inversion parameters 
the spherical coordinates Q = (pa, «a, Pb, «b, ^b) at two chosen observation 
dates A and B are used. Multivariate Gaussian proposal distribution Pt(Q',Qt) 
(describing the transition probability from last accepted spherical coordinates Q t 
to candidate spherical coordinates Q') is adopted for candidate orbit sampling. In 
particular, here we use multi-normal proposal distribution with highly correlated 
ranges pA and ps, as such a correlation is very often for short observa- 

tional arcs [15] (also see Fig. [I]). As the number of observations and observational 
time interval increases, the ranges distribution become more confined and less cor- 
related. This simple proposal distribution in spherical coordinates phase-space 
corresponds to highly complicated proposal distributions in Cartesian or/and Ke- 
plerian element phase space. To increase our chances of detection of multimodal 
distributions, we decided to use a number of shorter Markov chains instead of one 
long chain. We start our sampling with a set of random orbits (i.e., 10 sample 
orbits). Each of those composes a first element in a separate Markov chain. Next, 
a new candidate orbit is proposed (one candidate orbit for each chain) with a help 
of proposal densities. That orbit is then accepted into a chain or rejected according 
to MH acceptance criteria: 



If a> 1, then P 



t+i — 



P. 



If a< 1, 
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a 




(3) 
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J' and J t are the Jacobians from topocentric coordinates to orbital parameters 
for the candidate and the last accepted sample respectively. P' and P t denote a 
set of orbital elements, and p p is the a posteriori probability density function. In 
practice, this means that a trial orbit is always accepted if it produces a better fit 
to the full observational data set than the last accepted one in a given chain. The 
trial orbit is sometimes accepted if the fit is worse than the last accepted solution. 
If we have a choice between jumping to 'bad fit' and a 'very bad fit', the 'bad 
fit' jump is more probable. If a candidate orbit is rejected, then the last accepted 
solution is repeated in a chain. For more details on MCMC ranging please see |15j . 
To avoid manual tuning of the proposal density covariance, we have implemented 
an automated version of the MCMC algorithm. In the automated version, we 
first run 'pre-runs' (which could be considered as a burn-in period) to the actual 
sampling, in order to find adequate proposal density. That proposal density is 
computed based on empirical samples from the pre-runs, and then updated in a 
manner similar to that described in [8]: 

P?(Q' I Qi.-.Qn) ~ iV(Q t ,c^R« + c^I), (4) 

where R,W is the empirical 6x6 covariance matrix obtained from pre-run i from 
n samples Qi, Q ra ; c 2 d = 2.4/ y/d is a scaling factor, that depends on the number 
of dimensions d in the problem, in our case d — 6, Cd ~ 0.97. This scaling factor 
is optimized for for MH in the case of Gaussian targets and Gaussian proposals 
(For more details, please see Gelman et al. 1996 [1]). We stop updating the 
proposal covariance matrix after the pre-runs in order to avoid introducing biases 
in the sampling and to preserve the desired properties of Markov chains. The 
number of pre-runs depends on, for example, convergence diagnostic checks. The 
minimum number of pre-runs is two. The first proposal covariance matrix has 
to be user defined. Then, we start the automated MCMC ranging with the first 
pre-run and require 500 sample orbits. From that first pre-run, we compute the 
new proposal covariance matrix and continue with a second pre-run (1000 required 
orbits). A third pre-run is possible if the final number of required orbits is 50000 or 
more. After those pre-runs, we continue with iterative runs until all of the below 
conditions are fulfilled: 

• the algorithm has decided that it has reached the stationary distribution 
(please see the next section for more details), 

• the jump ratio in the last run is between 15-50% (recommended MCMC 
ratio), 

• the difference in the obtained maximum likelihood p.d.f in subsequent runs 
is smaller than a threshold value (| 2.0 log pdfi/pdfi^i |< 2.0). 
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• the difference in the obtained empirical covariance matrix in sub-sequentive 
runs is smaller than a threshold value (| 2.01ogi?j/i?j_i |< 2.0). 




0.(X>5 



Figure 1: Distribution of topocentric ranges pa and ps computed using a two-body model with 
six observations, 0.3-arcsec astrometric noise assumption, and Jeffreys' a priori. 

Even though the conditions are designed to ensure that we have reached the 
stationary distribution, one should always remember that convergence might still 
not have been reached (see the next section). If those conditions were not fulfilled 
and the number of iterations exceeded the maximum number of iterations (here 
we use five iterations after the pre-runs), the program exits with a warning that 
the stationary distribution has not been reached. Example results of a converged 
run are shown in Fig. |2j 
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Figure 2: Distribution of orbital elements (here: Cartesian ecliptic), computed using two-body 
model with six observations, 0.3-arcsec astrometric noise assumption and Jeffreys' a priori. 

2.1.2. MCMC ranging convergence diagnostics 

After obtaining sets of possible orbit solutions for each chain, we perform con- 
vergence diagnostics, which is a critical part of MCMC algorithms. Convergence 
diagnostics relate to the idea that a Markov chain, after a sufficient number of 
iterations, will eventually converge to the stationary distribution (target distribu- 
tion), starting from any point in the phase space and then mixing in that part 
of phase space forever. The proof of convergence can be found for example in 
Roberts and Rosenthal 2004 [17]. The main concerns are: (1) convergence to the 
stationary distribution; (2) the length of the burn-in period (for example, see Fig. 
[5]); and, finally, (3) number of iterations it will take to summarize the posterior 
distribution [2]. 

In general, there are two approaches to monitoring the convergence: monitoring 
of one chain for a long time or monitoring many chains for shorter periods of 
time. At present, a set of convergence diagnostics methods is available for MCMC 
algorithms for both of these approaches. The main idea usually involves running 
the Markov chain repeatedly from different initial states and checking if the chains 
all converge to approximately the same distribution or, in the case of a single 
chain, running it for long time and checking if it "looks stable". The diagnostics 
help identify problems with convergence, but cannot prove that the convergence 



6 




Figure 3: MCMC ranging with Jeffreys' a priori concentrates on the orbits with smalier Cartesian 
ecliptic coordinate x. 



has occurred. 

The basic convergence testing methods include history plots (plot of sampled 
parameters against the number of iterations, also called trace plots) and marginal 
distribution plots. If the algorithm converged to a stable distribution, the history 
plot will move around the mode of the distribution, and all the chains will cover 
the similar parameter range. Trending in trace plots can indicate lack of conver- 
gence. Figure [6] shows an example trace plot. If the separate chains converged to 
different areas in the phase space, we may be able also to see "lumpy" posterior 
marginal distributions, with distinct solution clouds. This effect can be sometimes 
also seen in the observed-minus-computed (O-C) residuals plots. To test the ran- 
domness of a data set, autocorrelation function plots (correlograms) are often used. 
Correlograms involve computing the correlation between data sets from a single 
chain and separated by some lag. If random, all the correlations will be near zero. 
Gelman-Rubin [5] proposed a statistical indicator, a so-called "shrink factor" R, 
involving the average within-chain variances, and the variance between the means 
obtained from m parallel over-dispersed chains. Once convergence is reached, the 
between chains and within chains variance should be almost equivalent because 
variation within the chains and variations between the chains should coincide, so 
R should approximately equal one. Rafery and Lewis [16J proposed a method for 
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Figure 4: MCMC ranging with uniform a priori concentrates on the orbits with larger Cartesian 
ecliptic coordinate x. 

estimating the total numbers of iterations needed, the length of the burn-in period 
based on chosen quantile, required accuracy, required probability to obtain speci- 
fied accuracy and convergence tolerance. Gewke [6] indicated the use of methods 
from spectral analysis in testing for convergence. The method is graphical and 
based on running means of parameters and can be used in, for example, checking 
for the length of the burn-in period. Zellner and Min [22] proposed a criterion 
based on conditional probability to test for convergence to a correct distribution. 
Extensive review of MCMC convergence diagnostics methods can be found in for 
example Mengersen et al. 1999 and Cowles and Carlin 1996 [EJE]. Topic of biases 
induced in the MCMC diagnostics is considered in for example Clowes et al. 1999 
[TTj . In our analysis, we make use of trace plots, kernel plots, 0-C residuals plots, 
and Gelman-Rubin diagnostics. 

2.2. Collision probability 

To assess the impact probabilities, we use techniques that rely on the orbital- 
element probability density function (p.d.f.) characterized using MCMC orbital 
ranging [T51IT] and MC ranging [20]. After obtaining a set of possible orbit solutions 
(here we use 5 x 10 4 different MCMC orbits for all the objects or/and 5 x 10 4 MC 
orbits), we propagate each of the orbits into the future to check for a possible 
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Figure 5: The marginal distribution of ecliptic coordinates i and x, computed for asteroid 2008 
TC3 using five observations and 0.5-arcsec astrometric noise assumption. The arrow indicates 
the burn-in samples, before reaching the stationary distribution. 
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Figure 6: A history plot generated from MCMC sampling of orbits for asteroid 2008 TC3, 
assuming 1.5-arcsec astrometric noise and including four observations for a time interval of 0.003 
day. 
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collision with the Earth within the determined impact interval r (time interval 
when collisions are possible). We step through this interval depending on the 
approximated time needed for the object to hit the Earth, given as it would be 
heading straightforward to the Earth. The next propagation step is therefore 
fixed as the time interval needed for the object to reach the Earth, given it would 
be heading head-on with our planet and depending on the current velocity and 
distance. The collision probability is given by [THl E] : 

P c (t) = [ dp p p (p,t), (5) 
Jo 

where R is the radius of the planetary body and p p (p, r) is defined as: 

Pp (p',r) = J dPp p (P)6 D (p' - p(P,r)), (6) 

where p(P, r) is the minimum distance of the small body from the center of the 
large (planetary) body, 5d is Dirac's function, P denotes orbital elements, p p is the 
a posteriori p.d.f. of the orbital elements leading to a collision. For more details, 
please see Virtanen and Muinonen 2006 and Muinonen et al. 2001 [191 EH]- I n 
case of 2008 TC3, we used an impact interval of Oct. 6-8, 2008. In the case of MC 
results, the orbits are weighted with the a posteriori p.d.f.; in the case of MCMC 
results, the weight of an orbit equals the number of repetitions of the orbit. 

3. Results and discussion 

To study the evolution of collision probabilities for the asteroid 2008 TC3 based 
on a short observation arc, we use a limited number of observations from the 
discovery night only. Here we use the first six observations obtained on October 
6, 2008. We divide this observational data into separate small sets. We start with 
two observations and then the subsequent sets contain one more observation than 
the previous one. We end up with five sets of different length (please see Table |l| 
for details). 



Number of Observations 


2 


3 


4 


5 


6 


Time interval [days] 


0.01 


0.02 


0.03 


0.04 


0.06 



Table 1: Observational data sets. 



First, we obtain a set of MCMC and MC orbits using each of those sets, and 
then we perform the collision probability computation. We repeat this process for 
different assumptions of the observational noise and different dynamical models. 
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Figure 7: Time evolution of collision probabilities for 2008 TC3 using MCMC ranging generated 
orbits with Jeffreys' a priori (n-body approach). 
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(a) Green - the marginal distribution of the orbital elements 
obtained for asteroid 2008 TC3 using five observations and 
0.5-arcsec noise. Blue - orbits leading to a collision with the 
Earth from the green set. Red - orbits leading to a collision 
with Earth obtained using 3 observations and 2.0-arcsec noise. 
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(b) Green - the marginal distribution of the orbital elements 
obtained for asteroid 2008 TC3 using 6 observations and 0.3 
arcsec noise assumption. Red - orbits leading to a collision 
with the Earth obtained using 3 observations and 2.0-arcsec 
noise. 



Figure 8: Shift to non-impact zone in the phase space in the two-body approximation (adding 
the 6th observation leads to a vanishing collision probability) . The black dot indicates the least- 
squares solution obtained with the full observational data set (see table [7]). This orbit does 
not lead to a collision with the Earth in the specified time interval when propagated using the 
two-body model. (For interpretation of the references to color in this figure caption, the reader 
is referred to the web version of this article.) 
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(a) Green - the marginal distribution of the orbital elements 
obtained for asteroid 2008 TC3 using five observations and 
0.5-arcsec noise. Blue - orbits leading to a collision with the 
Earth from the green set. Red - orbits leading to a collision 
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(b) Green - the marginal distribution of the orbital elements 
obtained for asteroid 2008 TC3 using 6 observations and 0.3- 
arcsec noise. Blue - orbits leading to a collision with the Earth 
from the green set. Red - orbits leading to a collision with the 
Earth obtained using 3 observations and 2.0-arcsec noise. 



Figure 9: Shift to impact phase space in the n-body approach. The black dot indicates the least- 
squares solution obtained with the full observational data set (see table [7]) . (For interpretation 
of the references to color in this figure caption, the reader is referred to the web version of this 
article.) 
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Nr. of 


Noise 


MC 


MCMC 


obs. 


[arcsec] : 


ranging: 


ranging: 






Jeffrey's apriori 


Jeffrey's apriori 




0.3 


0.75 (1536) 


0.58 (12587) 




0.5 


0.62 (1722) 


0.57 (12157) 


2 


1.0 


0.65 (2731) 


0.61 (10815) 




1.5 


0.67 (2586) 


0.61 (11136) 




2.0 


0.69 (2274) 


0.62 (11124) 




0.3 


0.008 (178) 


0.008 (308) 




0.5 


0.02 (282) 


0.01 (490) 


3 


1.0 


0.14 (663) 


0.11 (1218) 




1.5 


0.54 (1403) 


0.22 (3062) 




2.0 


0.39 (3190) 


0.28 (4521) 




0.3 


0.0007 (191) 


0.0006 (83) 




0.5 


0.005 (264) 


0.004 (308) 


4 


1.0 


0.02 (388) 


0.01 (725) 




1.5 


0.04 (591) 


0.02 (1042) 




2.0 


0.04 (967) 


0.06 (1666) 




0.3 


0(0) 


0(0) 




0.5 


0.0002 (325) 


0.0001 (11) 


5 


1.0 


0.007 (414) 


0.005 (457) 




1.5 


0.007 (516) 


0.01 (739) 




2.0 


0.02 (566) 


0.02 (1065) 




0.3 


0(0) 


0(0) 




0.5 


0.0000008 (2) 


0(0) 


6 


1.0 


0.0005 (449) 


0.0006 (63) 




1.5 


0.004 (490) 


0.004 (417) 




2.0 


0.01 (568) 


0.009 (694) 



Table 2: Collision probability, two-body approach. Number of different collision orbits (out of 
50 000 different orbit solutions) given in brackets. 

The results obtained are listed in Table [2] for the two-body case and in Table [3] for 
the n-body case. 

We list the collision probability computed for the two-body approximation for 
MC and MCMC ranging in Table [2] for comparison between the two methods. 
The results obtained from both methods agree very well. In Table |3j we list the 
impact probability computed in the n-body approximation, for different a priori 
p.d.f. assumed. The collision probability computed with the n-body and Jeffreys 
a priori results in a high collision probability already starting with 2 observations, 
and continues to grow with subsequent observations added reaching approximately 
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Nr. of 


Noise 


MCMC 


obs. 


[arcsec] : 


ranging: 






Jeffrey's apriori 




0.3 


0.80 (15281) 




0.5 


0.79 (15114) 


2 


1.0 


0.80 (14584) 




1.5 


0.79 (13991) 




2.0 


0.79 (13585) 




0.3 


0.61 (2813) 




0.5 


0.63 (3271) 


3 


1.0 


0.77 (6025) 




1.5 


0.80 (8412) 




2.0 


0.82 (12113) 




0.3 


0.51 (14960) 




0.5 


0.44 (7878) 


4 


1.0 


0.51 (4742) 




1.5 


0.63 (6116) 




2.0 


0.70 (8678) 




0.3 


0.92 (44053) 




0.5 


0.71 (28966) 


5 


1.0 


0.53 (12620) 




1.5 


0.53 (9611) 




2.0 


0.57 (8930) 




0.3 


« 1.0 (49771) 




0.5 


0.95 (46079) 


6 


1.0 


0.69 (27786) 




1.5 


0.59 (16348) 




2.0 


0.56 (11939) 



Tabic 3: Collision probability (as a weighted fraction of collision orbits), n-body approach. 
Number of different collision orbits (out of 50 000 different orbit solutions) given in brackets. 
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1.0 for 6 observations and 0.3 arcsec astrometric noise. Adding additional obser- 
vations (above 6 observations) to the data sets used leads to a collision probability 
approaching 1.0. Collision orbits are located in a small phase-space area, as com- 
pared to the phase space of all permitted orbital solutions, and have relatively 
high weight factors. The usage of Jeffreys' a priori shifts the solutions area to 
concentrate more on the orbits with smaller semi-major axes, whereas the usage 
of a uniform a priori causes a shift towards larger semi-major axes orbits (please 
see Figs. [3] and [4] for comparison of marginal distributions for orbital elements 
using Jeffreys' and uniform a priori). Results should be tested against different a 
priori p.d.f. assumptions. 

In the case of short observational arcs, the topocentric ranges at the observation 
dates A and B are very often highly correlated. An example of this can be found in 
Fig. [TJ which presents the distribution of topocentric ranges pa and ps computed 
using a two-body model, six observations, 0.3 arcsec astrometric noise, and Jeffreys 
a priori. The usage of a two-body dynamical model led to a shift of solutions to 
a non-collision phase space (see Fig. [8J, whereas the usage of the n-body model 
led to a shift into the collision phase space (see Fig. 9|. For comparison, we also 
computed a least-squares solution, please see Table [7 

Both methods show a sensitivity to the noise assumption. In general, the sen- 
sitivity to the noise should become small for extensive observational arcs (the data 
contains more information). However, in our case (maximum six observations and 
0.006-day time interval), such a stability cannot be yet observed. The sensitivity 
of the derived probabilities to the noise assumption can be estimated as [HI EE] : 



Sensitivity 



AlogPc 
Act 



(7) 



Large nonzero values indicate that the results are highly sensitive to the noise 
assumption. In Tables [4] and [5j we list the sensitivity estimation (per arc sec) for 
collision computation for all the sets, for both methods, using the two-body and 
n-body approaches. 



Obs. arc [days] 


MC 2-body 


MCMC 2-body 




Jeffreys apriori 


Jeffreys apriori 


0.01 


0.02 


0.02 


0.02 


1.8 


1.0 


0.03 


0.99 


0.7 


0.04 


1.33 


1.53 


0.06 


2.73 


1.18 



Table 4: Sensitivity to the noise assumption for MC and MCMC ranging techniques (2-body 
approach) . 
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Obs. arc [days] 


MMCC n-body 




Jeffreys apriori 


0.01 


0.0032 


0.02 


0.08 


0.03 


0.13 


0.04 


0.12 


0.06 


0.15 



Table 5: Sensitivity to the noise assumption for MC and MCMC ranging techniques (n-body 
approach) . 



Nr. of 


MCMC 


obs. 


ranging: 


2 


0(0) 


3 


0.43 (43491) 


4 


0.70 (55407) 


5 


0.55 (35895) 


6 


0.54 (35829) 



Table 6: Collision probability (as a weighted fraction of collision orbits) computed with astro- 
metric noise sampling, n-body dynamical model, and Jeffrey's apriori used. Number of different 
collision orbits (out of 50 000 different orbit solutions) given in brackets. 

We indicate the need for more reliable methods of collision probability esti- 
mation (less sensitive to the noise assumption) for poorly observed objects. Here 
we also perform collision probability computation using MCMC computed orbits, 
where astrometric noise has been included in the list of sampling parameters. Ta- 
ble [6] shows a result of that exercise. The time evolution of this probability differs 
from the one computed using a fixed assumed astrometric error. The impact prob- 
ability also starts with a quite high value, then it grows, drops down, stabilizes 
and grows again (for data sets exceeding six observations). We have assumed that 
all the observations have the same astrometric uncertainties. All the observations 
that we used come from a single survey, that is the Mt. Lemmon Survey (1.5-m 
telescope). We have selected a priori information on the astrometric uncertainties 
at Mt. Lemmon based on astrometric residuals of all objects observed in 2008 
at that observatory. We have assumed a Gaussian a priori for the noise with a 
covariance matrix computed from those residuals, and mean of 0.3 arc seconds. 
One could go even further by sampling the astrometric accuracy for each observa- 
tion separately, for each observatory or for each observer, and consequently using 
different a priori information for those. This kind of analysis is, however, beyond 
the scope of this paper. 



18 



Element 


Value 


1-cr variation 


Units 


X 


0.978354962 


0.57E-07 


AU 


y 


0.2242293386 


0.45E-08 


AU 


z 


0.871659598E-03 


0.11E-07 


AU 


X 


-0.776631371E-02 


0.55E-07 


AU/day 


y 


0.1720023476E-01 


0.45E-08 


AU/day 


z 


-0.755199990E-03 


0.11E-07 


AU/day 



Table 7: Least square solution (cartesian ecliptic coordinates) computed with the complete set 
of observations 859 (70 discarded) for the epoch of 2008/10/6.0 using OpenOrb software [7J. 1.0 
arcsec astrometric noise assumption and n-body dynamical model. 

4. Conclusions 

We have computed and examined collision probabilities using the MC and 
MCMC ranging methods for the Earth-impactor asteroid 2008 TC3 using small 
amounts of data. We have considered sampling of the astrometric uncertainty as 
an additional parameter to the orbital elements, which can be important for cases 
where that uncertainty is not known or cannot be assumed. We indicate the need 
for more reliable methods of collision probability estimation (less sensitive to the 
noise assumption) for poorly observed objects. New ground-based and satellite 
sky surveys (like the forthcoming Gaia survey), due to their improved accuracy, 
will significantly improve the computation of collision probabilities based on small 
amounts of data. NEO impact monitoring is currently not included in the Gaia 
processing pipeline and will have to be done outside the main data processing. 
Fast access to the Gaia data is essential for objects such as the Earth impactor 
2008 TC3. Earlier data realizes for NEOs should be therefore considered. 
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